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We study the properties of discrete breathers, also known 
as intrinsic localized modes, in the one-dimensional Frenkel- 
Kontorova lattice of oscillators subject to damping and ex- 
ternal force. The system is studied in the whole range of 
values of the coupling parameter, from C = (uncoupled 
limit) up to values close to the continuum limit (forced and 
damped sine-Gordon model). As this parameter is varied, the 
existence of different bifurcations is investigated numerically. 
Using Floquet spectral analysis, we give a complete character- 
ization of the most relevant bifurcations, and we find (spatial) 
symmetry-breaking bifurcations which are linked to breather 
mobility, just as it was found in Hamiltonian systems by other 
authors. In this way moving breathers are shown to exist even 
at remarkably high levels of discreteness. We study mobile 
breathers and characterize them in terms of the phonon radi- 
ation they emit, which explains successfully the way in which 
they interact. For instance, it is possible to form "bound 
states" of moving breathers, through the interaction of their 
phonon tails. Over all, both stationary and moving breathers 
are found to be generic localized states over large values of 
C, and they are shown to be robust against low temperature 
fluctuations. 

FAGS numbers: 05.45. -a,73.20.Ry 



I. INTRODUCTION 

The phenomenon of (non-topological) locahzation in 
discrete nonlinear lattices {i.e. intrinsic localized modes 
or discrete breathers) has received a great deal of at- 
tention from both theoretical and (as of lately) experi- 
mental research. Indeed, recent observations of dis- 
crete roto-breathers in Josephson-junction ladder circuits 
have placed the subject on a firm experimental footing 
(see also Q]). Most of the theoretical and computational 
work on discrete breathers has dealt with Hamiltonian 
systems, of fundamental interest in Physics. For a re- 
view, see Refs. (5]J^. Comparatively, the easier case of 
dissipative breathers has received much less attention, 
though the experimental systems that we have just men- 
tioned belong to this class. 

Mathematical proofs of existence of discrete breathers 
in rather general dissipative networks of oscillators piq] 
appeared soon after those of Hamiltonian networks | B| . 
While in the later condition of non-resonance of 



the localized oscillation with the band of extended nor- 
mal modes of the lattice has to be satisfied, that is not 
an issue in the forced-damped case, and the dissipative 
breather possesses the character of attractor for initial 
conditions in the corresponding basin of attraction. As 
archetypical example of Klein-Gordon lattices of oscilla- 
tors, we consider the standard Frenkel-Kontorova model 
with commensurability one {i.e. average interparticle dis- 
tance equal to the period of the sinusoidal substrate po- 
tential). In section || we discuss the numerical proce- 
dures used to obtain accurate breather solutions, which 
are based on the continuation from the uncoupled limit 
of the model. 



In section HI we explain some general features of the 



linear stability (Floquet) analysis of forced-damped peri- 
odic discrete breathers. For the sake of readability, this 
section is intended to be self-contained, to some extent. 
After deriving some straightforward properties of the Flo- 
quet multipliers, we obtain some formal conditions for the 
non-appearance of extended instabilities of the uniformly 
oscillating background, along with the tail analysis valid 
for not-too-large forcing. The good fitting of our numer- 
ical data to the results of this section ensures its validity 
for the parameters used in the numerical work. 

Section IV reports on our numerical findings, con- 



cerning pinned discrete breathers, which are summa- 
rized in the phase diagram against the coupling param- 
eter. Pitchfork and Andronov-Hopf bifurcations sepa- 
rating different periodic and quasiperiodic breathers ap- 
pear as generic features of this phase diagram. At very 
high values of the coupling parameter, when the width 
of the discrete breather is much larger than the period of 
the substrate potential, a Goldstone mode in the Floquet 
spectrum signals the approach to the continuum limit. 

In section ^ we study the mobility of discrete 
breathers, a subject which is yet poorly understood. Af- 
ter discussing the procedures used to obtain and continue 
mobile breathers we explain successfully, with the aid of 
simple physical arguments, the numerical power spectra 
in the tails. Then we study collisions between discrete 
breathers. We find all possible scenarios, ranging from 
"elastic" to completely "inelastic" collisions; this latter 
case includes both breather annihilation and, more in- 
terestingly, the formation of "breather molecules" which 
can be either pinned or mobile. Finally, in section VI, we 
summarize the main conclusions of our work. 
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II. MODEL AND BREATHER GENERATION 

The equations of motion of the Frenkel-Kontorova 
chain subject to damping and an (spatially uniform) ex- 
ternal driving force are, in dimensionless form, 

I . , 

ilj + aiij + — sin(27rUo) = 

C (uj+i - 2uj + Uj-i) + FacSin{ujbt) (1) 

In order to generate a discrete breather configuration 
we start in the anti-integrable (uncoupled) limit C = 
0, using two different amplitude attractors of the single 
pendulum equation of motion. That is, we first consider 
the dynamics of a single forced and damped pendulum, 
and try to find a region of parameters where there is 
at least two different attractors coexisting. Note that, 
generically, all oscillators have at least two attractors for 
sufficiently low values of the damping a and the force Fac, 
if the frequency of the force ujb is not wildly different from 
the typical frequencies of the autonomous oscillator. 

Therefore we initially choose values for a, Fac, and LUf,, 
and keep them fixed while we vary C. Then, for instance, 
we fix one of the oscillators to the high amplitude solu- 
tion and all the others to the low one. Using as initial 
condition this anti-integrable configuration, we turn on 
adiabatically the coupling parameter C. Following Sepul- 
chre and MacKay's work on dissipative breathers [Q, the 
initial solution can be continued for C 7^ 0, at least until 
a bifurcation is reached. That paper shows how, in con- 
trast to Hamiltonian systems, forced and damped sys- 
tems have it easier to comply with the conditions of the 
continuation theorem, since there is no problem of reso- 
nance with phonons (we have attractors), and the relative 
phases of the oscillators are locked by the external force. 

Moreover, if the variation in C is small enough, the 
discrete breather remains an attractor of the dynamics 
(since one expects the basins of attraction to evolve con- 
tinuously with C as well). This makes the numerical 
continuation greatly simpler: it is possible to just vary 
adiabatically the coupling C as we integrate the equa- 
tions of motion (|l|), and the dissipative dynamics drives 
the system to the stable attractor. Contrast this with 
the expensive root-finding methods one needs to use for 
breather continuation in Hamiltonian problems Q . 

In addition, we performed a linear stability analysis of 
the periodic solutions (Floquet-Bloch analysis, see below) 
in order to investigate the nature of possible bifurcations. 
In some cases we have added to the initial conditions a 
small random noise (typically of order 10~^) to test for 
robustness. We have taken special care dealing with fi- 
nite size effects. For low values of C (C < 0.6) small 
lattice sizes can be used (say N — 40). However, once 
the breather is dressed by a phonon tail (see below), we 
have needed to increase the lattice size up to = 900 
in order to avoid finite size effects. We think this is an 
important point to check since experiments in real dissi- 
pative systems are done in small lattices . Numerical 



integration of equations of motion has been done using 
a fourth-order Runge-Kutta scheme. Most of the simu- 
lations in this paper have been done with the following 
parameters: a = 0.02, ujf, = 0.27r, and Fac = 0.02, al- 
though we have sometimes changed them to confirm the 
general validity of our results. 

III. LINEAR STABILITY ANALYSIS 

A. Floquet multipliers 

Let us consider a small perturbation, {vj{t)}, of the 
discrete breather {uj{t)} solution, Vj — Uj -I- e^. After 
direct substitution in the equations of motion and dis- 
carding terms which are nonlinear in e^, one finds 

'ij + aij + cos(27rUj(<))ej = C (cj+i — 2ej + ej-i) (2) 

These form a system of coupled linear differential equa- 
tions with time periodic coefficients, for Uj{t) is a periodic 
function of time. For a system of size N , the integration 
of the linearized equations (|^) over a period tf, = 27r/wb 
of each of the 2N vectors {ej(0), £^(0)} forming some ba- 
sis of the tangent space defines the 2A'' x 2N Floquet (or 
monodromy) matrix J- 

that relates the small perturbations at t — ti, to those at 
t = 0; in other words, J- is the matrix associated to the 
tb-map of (||). 

The linear stability of the breather solution {uj{t)} 
requires that all the eigenvalues of the Floquet matrix 
(called also Floquet multipliers) are inside the unit circle. 
Since is real, if 11 is an eigenvalue of its complex 
conjugate Jl is also an eigenvalue of T. But the Floquet 
spectrum has more structure, since one can transform 
the linear system of Eq. (H) into a Hamiltonian one (see 
Rcf. [|o|). By transforming the ej variables according to 

e,it) = e-'^'/\{t), (4) 

this yields 

fjj — ^a^/4 — cos(27rMj(t))^?7j = 

Cfe+i - 2r?j +?7j_i) (5) 

These are the equations of motion of a (non autonomous) 
Hamiltonian system of oscillators, for which the eigenval- 
ues of the (symplectic) i^-map must come in pairs such 
that their product is unity. Together with the fact that 
the map is real, one has these well known three pos- 
sible cases: (i) pairs of eigenvalues lying on the unit cir- 
cle, with Ai = A2; (a) pairs lying on the real axis, with 
Ai = I/A2; (Hi) 4-tuples of eigenvalues with Ai = A3, 
A2 = A4, Ai = I/A4. 
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Since the transformation scales the eigenvalues by 
a factor exp[— Q!i;,/2], the Floquet multipliers of (§) must 
either lie on a circle of radius exp(— atb/2), or on the real 
axis such that ^ifi2 — exp(— at{,), or come as 4-tuples 
such that ^1 = /i3, /i2 = M4, Ail = exp{—atb) /jl^. 

An important difference with the Hamiltonian case, 
where a "phase" and the "growth" modes |^ are always 
associated to the double eigenvalue +1 in the Floquet 
matrix of the discrete breather, is that these modes do 
not exist for the forced-damped case. The reason for that 
is that both the breather frequency and the time origin 
are fixed by the external force, so that the associated 
degeneracies are removed. 

B. Extended instabilities 

In the limit of an infinite system (A'' — > oo), the spec- 
trum of J- consists of a continuous part associated with 
spatially extended eigenvectors and a discrete part asso- 
ciated with spatially localized eigenvectors. The contin- 
uous part of the spectrum of J- is the continuous spec- 
trum of the linearized problem around the homogeneous 
solution {i.e., without breather) Uuit)} = {uoo(t)}. As 
pointed out by Marfn and Aubry 02] , using the fact that 
the limit (in the appropriate sense) of the sequence of 
spatial translations of the Floquet matrix ^ of the sys- 
tem with breather is the Floquet matrix !Fo of the system 
without the breather, one proves easily that the spectrum 
of J-Q is included in the spectrum of J-. Reciprocally, 
the limit of the sequence of spatial translations of an ex- 
tended eigenvector of !F can be seen to belong to the 
spectrum of J-q. 

First, we are going to consider the spectrum of Tq, 
so we now will pay attention to the linearized equa- 
tion of motion around the homogeneous solution 
of (|l|), {uj{t)} — {uoo{t)} and denote simply f{t) = 
cos(27ruoo(^))- Under the usual periodic boundary con- 
ditions, we look for solutions of the linear problem with 
the plane- wave form 

Vjit) = e">^x''{t) (6) 

In other words, x''(0 the (spatial) Fourier coefficient 
of rjj[t). Inserting (^ into the equations (^), and denot- 
ing by E{q) — 4Csin^((7/2) — we have, for each 
value of (7, the equation 

X«(i) + (i?(g) + /(t))x^(<)-0 (7) 

This is a Hill equation. For each solution x''(t) of the 
single Hill equation (^ we have a solution of the form 
(^ for the equations (P), and thus, a solution 

for the linearized problem. The Hill equation (Q) has a 
general solution which can be expressed in terms of its 
normal solutions, which have the property 



x'(i + 27r/c.fc) = A,x«(t), (9) 

where is called the characteristic number of the 
equation. The complex number pq defined as Xq = 
exp{2'Kpq/u}b) is the called characteristic exponent (its 
imaginary part being defined up to an additive multiple 
of Wfc). In the generic case in which equation (Q) has two 
different characteristic numbers A^, A~, their product is 
equal to unity, A^j^A^ = 1, and the general solution has 
the form 

X'(i) - c+ef^'^^it) + C-eP^'^qit) (10) 

where c+,c_ are constants and 'ipq,^q are time peri- 
odic functions with period 27r/a;{,. Consequently, x'^{^) 
is bounded by A'exp(p™^'^t), with K some constant, 
and p^'^^ = max{p+,p^}. Thus, from equation (||), we 
conclude that the stability of the homogeneous solution 
{moo(0} is assured in the parameter region in which 

p,up = supp:;-'^<a/2 (11) 

The determination of this region in parameter space 
can only be made by numerical means. For the range of 
parameters that we have used in our study of damped- 
forced breathers, the function f{t) is a low amplitude 
oscillation around the value 1, and, as expected from 
the well-known results on weakly time dependent Hill 
equations, we have not observed instabilities by extended 
modes. 

In the next section, we will follow the continuation of 
the breather solution from the uncoupled limit, for in- 
creasing coupling and numerically compute the eigenval- 
ues of the Floquet matrix T. This will allow the charac- 
terization of the different bifurcations that the breather 
experiences when the coupling parameter increases. 

C. Tail analysis 

To proceed a bit further, we will assume from now 
on in this section that Mqo {t) is an oscillation of very low 
amplitude, so that for |j| ^1 the coefficient cos(27ruj(t)) 
in equations (||) is essentially unity, if one discards terms 
less than or equal to u^it). Then we are left with the 
standard problem of a linear chain with damping, which 
we can solve exactly (a similar analysis to the one below 
appears in Rcf. Q). 

Let us consider a semi-infinite chain with the boundary 
condition at the beginning given by eo(0 — exp(— iwt), 
and look for solutions of (g) of the form 

ej(t) = e(-«+*«)^'e-*'^* (12) 

Inserting ( |l2|) into equations (^) one obtains for the 
real and imaginary part, respectively 
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cosh^ cosq = 1 + 7-^(1 — ijJ^ 



sinh ^ sin q 



2C 



(13a) 
(13b) 



First, we will analyze the situation in which a = 0, and 
see how one recovers the well-known results for Hamilto- 
nian discrete breathers 15|. For ^ = 0, one has the 
familiar normal mode solutions, where the frequencies are 
given by the dispersion relation 



= 1 + 4Csin2(q/2) 



(14) 



and q (— tt < q < tt) is the wave-vector of the nor- 
mal mode. It is customary to denote loosely the nor- 
mal modes as "phonons" , and the interval of values of uj 
defined by (|l^ as "phonon band" . For ^ ^ one has 
exponentially decaying solutions 



,{t) = (-l)^e-«^e-^"* 



(15a) 
(15b) 



where the inverse decay length ^ and the frequency oj are 
related, respectively, through: 



= 1 -4Csinh2(^/2) 
= 1 + 4Ccosh2(^/2) 



(16a) 
(16b) 



Note that the values of lo in (^6|) are, respectively, 
below and above the phonon band, so we observe how 
the Hamiltonian linear lattice damps out any solution 
with a frequency component outside the phonon band 
(p^, while the normal modes are extended (^ = 0). As 
a consequence, a Hamiltonian breather needs to have 
all breather harmonics nuji, out of the phonon band, 
and then they decay exponentially with the character- 
istic length £_~^{nujb)- Thus the size of the Hamiltonian 
breather is given by = sup„ ^~^(nwb). 

When a 7^ 0, we have that ^(w) 7^ 0. Thus, any so- 
lution decays exponentially. For very low values of the 
damping, and frequencies well inside the (Hamiltonian) 
phonon band, the deca y len gth f is very large, so that 
sinh^ ~ ^ in equation ( |l3b|) , and thus it can be approxi- 
mated by 



2Csmq 2vg 

auj a 



(17) 



where Vg = (duj/dq) is the "group velocity" of the cor- 
responding normal mode, obtained from the dispersion 
relation (H). This approximation admits a simple phys- 
ical interpretation in terms of the competition between 
the damping and the velocity Vg at which the wave gen- 
erated (at the beginning of the semi-infinite chain) by 
the sustained perturbation propagates: The amplitude 
of the excited phonon decays in time as exp(— at/2), so 
that the time after which the amplitude has decayed by 
a factor of (1/e) is 2/ a, and thus the distance traveled 
by the phonon is 2vg/a. 



An example of the solutions ^(w) and q{uj) of equa- 
tions ([1^), for the particular values a = 0.02 and C = 
0.75, appears in Fig. |^. For comparison purposes, the 
graphs corresponding to the same value of coupling for 
the Hamiltonian case are included. 

Note that for the existence of damped-forced discrete 
breathers there is no need of a non-resonance condition 
(in contrast with the Hamiltonian case), because for any 
frequency lo, ^(uj) / 0. However, for low values of a, if 
some breather harmonic nwb belong to the interval of val- 
ues of UJ for which ^(w) is very small, the breather profile 
will show large "wings". In Fig. ^ we plot ^(Swb) as a 
function of the coupling C, for a = 0.02 and u)b = 0.2tt. 
Observe the dramatic decay of C(3ijJ6) at around C — 0.6 
corresponding to the entrance of the third breather har- 
monic in the (so to speak) "phonon band" , and compare 
the breather profiles for two values of C, respectively be- 
low and above, in Fig. |[ Both the wave-vector and the 
size of the wings in figure fit very well with g(3wb) and 
^(Swfc) from equations (O). 



IV. BIFURCATIONS AND PHASE DIAGRAM 

We have continued numerically breather solutions from 
the uncoupled (or anti-integrable) limit, for fixed values 
of the damping coefficient (a = 0.02), external force fre- 
quency {ujb = 0.27r) and intensity {Fac = 0.02). The spec- 
trum of Floquet multipliers was also numerically com- 
puted for each solution, thus monitoring their evolution 
on the complex plane. The configurations we have fo- 
cused on are two: the one-site breather and the two-site 
breather (adjacent sites). Of course many other config- 
urations are possible, by choosing from all combinations 
of sites in either the high-amplitude or low-amplitude at- 
tractor. However these two simplest breathers already 
provide a quite rich behavior, and, surprisingly, they al- 
low continuation into very high values of C, where the 
continuum limit is approached. 

At very low values of C the continued one-site breather 
is very narrow, symmetric around its localization site {i.e. 
U-i{t) — Ui{t) for all t and i), and all the Floquet multi- 
pliers lie on the circle of radius exp(— aifc/2) in the com- 
plex plane. The breather remains stable for increasing 
couphng up to the value Cpi = 0.52962 where a Floquet 
multiplier, which had previously detached along the real 
axis from that inner circle, reaches the unit circle at -1-1. 
The corresponding eigenvector of the Floquet matrix is 
localized around the breather site and possesses odd mir- 
ror symmetry with respect to that site, as shown in Fig. ||. 
Past the bifurcation, we are left with an unstable sym- 
metric breather and two new stable breathers, spatially 
asymmetric and one being the mirror image of the other. 
We can conclude that this is a forward pitchfork bifur- 
cation ]l6[ |, associated to a spatial symmetry-breaking 
transition of the discrete breather. 

In order to visualize the mirror symmetry breaking 
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character of this bifurcation, we plot in Fig. ^ the dif- 
ference A(0) = u_i(0) — Mi(0) at time t = (mod tt) 
of the positions of the neighbor oscillators on both sides 
of the localization site, as a function of the coupling 
parameter C in the vicinity of the bifurcation value. 
It is not surprising that, close to the bifurcation, the 
difference A(0) scales with the coupling parameter as 
(C — Cpi)^/^, because the one-dimensional character of 
the unstable manifold of the symmetric breather allows 
to reduce the analysis to that of a pitchfork bifurcation 
in a one-dimensional map, where this scaling behavior 
for the distance between branches is well-known. 

The asymmetric stable branches born at Cpi can be 
continued for higher values of coupling. We then observe 
how the amplitude of one of the neighbors of the central 
site keeps increasing, until it equals the amplitude of the 
central oscillator (which has in turn decreased slightly), 
at Cp2 = 0.55315. It should be noted that the rela- 
tive phases of the oscillators do not appear to change. 
At this value, a backward pitchfork bifurcation occurs, 
where the two stable asymmetric breathers and one un- 
stable symmetric two-site breather merge, and a stable 
mirror-symmetric two-site breather comes out. As in the 
bifurcation analyzed before, only one Floquet eigenvec- 
tor, associated to a Floquet multiplier of value -1-1, is 
involved. 

This unstable two-site breather which joins this sec- 
ond pitchfork is nothing but the two-site breather which 
can be continued from the uncoupled limit. Thus we have 
that the two elementary breathers constructed at the un- 
coupled limit, one-site and two-site, undergo an exchange 
of stability via this symmetry-breaking pitchfork mech- 
anism. For C < Cpi the one-site breather is stable and 
the two-site one unstable. For Cpi < C < Cp2 both 
are unstable, and the new asymmetric breather is stable. 
Past Cp2, the two-site breather is stable and the one-site 
unstable. This is exactly the same mechanism that was 
previously found by other authors [T^ for Hamiltonian 
breathers, and which is related to breather mobility as 
we explore in the next section. It is therefore plausible 
to conjecture that the mechanism is highly generic and 
might be expected in a large class of models. 

To be thorough in our description, a much less inter- 
esting bifurcation does appear in the two-site breather 
branch at very low C. The two-site breather is ini- 
tially stable at the uncoupled limit, but loses stability 
after a pitchfork bifurcation at C « 0.02. This is also 
a symmetry-breaking bifurcation as above, however it is 
instructive to investigate the differences: this time the 
spatial symmetry is broken in such a way that the new 
stable, asymmetric breathers suffer a de-phasing between 
the two central sites. This can be confirmed rigorously 
by examination of the relevant eigenvector at the bifur- 
cation. If one takes as reference for the time origin the 
instant at which the two central sites have maximum am- 
plitude, the unstable eigenvector for this second pitchfork 
shows components only in the velocity part, not in am- 
plitudes; the former bifurcation shows exactly the op- 



posite behavior. In any case, the continuation of the 
asymmetric breather from this bifurcation at C = 0.02 
is quickly lost after a Hopf bifurcation, and we have not 
found any more interesting behavior arising from these 
curious branches. 

Now we turn on to the continuation of the stable two- 
site breather branch past C > Cp2- We find now that 
near C = 0.817 two complex conjugate Floquet multipli- 
ers cross the unit circle at exp(±i(/3), with (p = 1.285, and 
the (periodic) two-site breather becomes unstable. The 
real and imaginary components of the associated eigen- 
vectors are shown in Fig. ^. Close to the bifurcation, 
small perturbations of the unstable periodic breather 
bring it into a quasiperiodic breather (as verified by in- 
spection of the Poincare section), thus confirming the 
scenario of an (Andronov-) Hopf bifurcation psf . In- 
deed, the power spectrum analysis of the quasiperiodic 
attractor (see figure) reveals two basic frequencies, Ub 
and uj-cicw — 0.0873. However, since this new frequency is 
rather different from the frequency ipu>b/2'K = 0.1285 as- 
sociated to the destabilizing eigenvalue couple, one con- 
cludes that the stable quasiperiodic attractor does not 
come out straight from the bifurcation. In other words, 
the simplest scenario of a subcritical Hopf bifurcation is 
discarded. Moreover, the quasiperiodic attractor can be 
easily continued back for lower values of the coupling 
{i.e., C < 0.817), therefore confirming that we have a 
subcritical Hopf bifurcation. Note that, as parameters 
other than C change, it is possible for these subcritical 
Hopf bifurcations to become supercritical, for their gener- 
icity lies in the Hopf character, not in being subcritical 
or supercritical. 

The lack of periodicity for the new (quasiperiodic) 
breather attractor prevents the use of Floquet analy- 
sis. However, its stable character can be numerically 
ascertained, checking its robustness against small per- 
turbations in the dynamics. This quasiperiodic two- 
site breather turns out to be stable for couplings up to 
C = 0.88, beyond which it starts moving spontaneously. 
Only when we reach C = 0.96 we recover again a stable, 
pinned quasiperiodic breather. Meanwhile, the pinned, 
periodic two-site breather (which became unstable after 
the Hopf bifurcation at C = 0.817) can be continued by 
a Newton method. At a value near C « 0.995, it rejoins 
the quasiperiodic two-site breather in an inverse (now 
supercritical) Hopf bifurcation, becoming stable again. 

We have tried to summarize most of these bifurcations 
in the sketch of Fig. H. We postpone to the next sec- 
tion the analysis of the observed mobile breathers in this 
region (and others) of parameter space. 

To conclude this section we will comment on the con- 
tinuation and stability analysis as C — > oo, i.e. the so- 
called continuum limit. The first interesting fact is that 
both the one-site and two-site breather, whether stable or 
unstable, have been found to be continuable for couplings 
as high as desired. In other words, they never disappear 
by, say, saddle-node bifurcations or the like. Another 
interesting point is that for C > 1 we have also found 
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pitchforks which connect the two branches, exactly in 
the same way as the first one at Cpi, Cp2- The ranges of 
couphng between forward and backward pitchfork bifur- 
cations (that is, where the connecting branches of stable 
asymmetric breathers exist) get progressively narrower 
for higher coupling values. And finally, both the breather 
profiles and their Floquet spectrum reveal a very natu- 
ral approach to the continuum limit: the solutions get 
broader in size, while the eigenvalue responsible for the 
pitchforks remains closer and closer to -1-1 at all times, 
announcing the appearance of the Goldstone mode (due 
to translational invariance) as C oo. These effects 
were ostensibly manifest at C ^ 5. 

V. MOBILE BREATHERS 

The problem of the mobility of discrete breathers is still 
very poorly understood. While the fundamental theory 
of stationary breathers is now firmly established |^ , mov- 
ing discrete breathers (MB for short) have so far eluded a 
rigorous treatment. But the fact is that moving breathers 
have been observed and studied through numerical sim- 
ulation in various works [^^|~ ^ , pT|j2^ ] , and they appear 
to be a phenomenon with the same degree of genericity 
as the stationary case. Up to now most of those works 
have dealt with Hamiltonian systems; here we present a 
study of moving breathers in our forced and damped F-K 
model. Some of the results are strikingly similar to those 
observed in Hamiltonian systems. 

We should first point out that a moving breather is 
something to be distinguished from a similar class of so- 
lutions, namely lattice solitons [ p3[|2^ ]. In a lattice soli- 
ton a pulse propagates without dispersion through the 
lattice, but, unlike breathers, there is no "internal" oscil- 
lation. This additional degree of freedom makes the mov- 
ing breather a more complicated object. For instance, in 
Hamiltonian lattices, it is easy to see that inevitably the 
moving breather resonates with the phonon band (note 
the presence of a quasi-periodic spectrum due to the addi- 
tional frequency introduced by the translational motion), 
and therefore it is not possible to have tails which decay 
to zero. It is not clear whether the solution is just a 
transient which eventually decays by phonon radiation, 
or maybe an infinite lifetime breather which "rides" on 
an infinite, small amplitude radiation background. 

But our model here is dissipative and has external forc- 
ing, and it turns out that moving breathers appear as 
proper attractors of the dynamics. Since these solutions 
are not transients, we can study and characterize them 
accurately and with great confidence. Even though this 
will not shed any light on the problem of existence of 
Hamiltonian moving breathers, there is another aspect 
of theoretical interest in which this study can contribute: 
the (possible) concept of a Peierls-Nabarro barrier for 
breather motion. This concept arises because of the sim- 
ilarities with the problem of mobility of discommensura- 



tions (kinks) in the Frenkel-Kontorova model . The 
discommensuration is an equilibrium static structure, for 
which one may ask how much energy it costs to displace 
it by one lattice site, until it reaches the equivalent con- 
figuration by (discrete) translational invariance. This is 
commonly referred to as the Peierls-Nabarro (PN) bar- 
rier. And it is possible to give a precise definition: from 
all possible continuous deformations of the initial con- 
figuration into the final one, take the one in which the 
maximal energy change along the path is the minimum. 
Very fruitful results in the theory of the F-K model have 
stemmed from this definition ]25t . 

But a corresponding definition of a PN barrier for 
breathers proves quite problematic. The difhculty lies in 
that it is not clear which space to use, since the configura- 
tions are now periodic functions, not static points. Some 
authors have suggested possible candidates for a rigorous 
definition, but the issue is still under debate |p^,p|,p7t . 
Technicalities aside, it is still possible to give a working 
definition of the Peierls-Nabarro barrier for breathers, 
at least in some cases. Most studies generate moving 
breathers by perturbing stationary ones, and Ref. p^ ] 
gave a systematic method to do this. Looking at the lin- 
ear stability analysis of the breather (Floquet analysis), 
they found that in many cases one can identify an eigen- 
mode which is distinctively localized and whose spatial 
symmetry is the appropriate for breather motion (note 
that similar depinning modes are responsible for the de- 
pinning of discommensurations ]2^ ] under uniform forc- 
ing) . It was found that adding a perturbation along this 
depinning eigenmode, provided one overcomes a certain 
threshold, results in a moving breather. Such thresholds 
are probably the best pragmatic approach to the defini- 
tion of a PN barrier. 

A further study p7| showed that many Hamiltonian 
lattices exhibit a very interesting behavior which is linked 
to mobility. It was found that, as the coupling is in- 
creased from C = 0, the one-site breather and the 
two-site breather each undergo a pitchfork bifurcation, 
where new branches of periodic but spatially asymmet- 
ric breathers emerge. These branches do in fact con- 
nect those pitchfork points, and the corresponding Flo- 
quet eigenmodes responsible for the bifurcations obvi- 
ously show a spatial symmetry which we could dub as 
"depinning" . This is exactly what we have found in our 
dissipative model, as shown in the previous section. And, 
just as in the cited works on Hamiltonian systems, we 
have also verified that the mobility is greatly enhanced 
for values of coupling in the vicinity of these bifurca- 
tions: very small amounts of perturbation along the de- 
pinning mode are enough to create the mobile breather. 
The upshot is that this phenomenon provides a mecha- 
nism for the existence of mobile breathers at relatively 
low couplings (high discreteness) and with very slow ve- 
locities, two properties which were counterintuitive and 
unexpected. 

We should note that the continuous sine-Gordon equa- 
tion under external ac forcing and losses does not sup- 
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port MB solutions The way in which a continuous 
breather destabiUzes is by a transition to a quasiperiodic 
state and finally creation of a kink-antikink pair | |30[ |. 

In the following we begin exploring this relation be- 
tween the stability of stationary breathers and the exis- 
tence of their mobile counterparts. Then we concentrate 
on studying the properties of moving breathers in dissi- 
pative systems in more detail. Finally, we explore other 
aspects such as collisions. 



A. Generation and phase diagram of MB 

We have found MB as proper attractors of the dy- 
namics in a wide range of couplings, in particular for 
0.5 ^ C < 0.96. We have found them either by excita- 
tion of the depinning mode of stationary breathers (as 
explained below) or simply by letting the system evolve 
to a steady state after the instabilities of some stationary 
breathers develop fully. Then it is possible to carry out 
a continuation of the MB into other parameter values, 
since their attractor property allows to change slowly a 
parameter and track the MB solution. Figure ^ shows 
the phase diagram of MB we have constructed with this 
procedure. 

Around the first symmetry-breaking bifurcation there 
is a narrow region in which we found MBs with regular 
motion and well defined velocity (Fig. A similar 

region is found in the interval 0.7 ;$ C < 0.88. We call 
these solutions induced fast breathers. To generate these 
steady state MB we have followed the procedure described 
in [ p8| . Surprisingly, this method works very well not 
only near the bifurcations. We typically use as initial 
conditions: 



u,{t = 0) = u° + Xef, 



(18) 



where is a stationary breather solution, ef corresponds 
to the antisymmetric (and localized) eigenvector mode at 
the bifurcation taking place near C = 0.53, and finally 
A measures the strength of the perturbation applied. It 
is found that, as in the Hamiltonian case, a critical Ac 
is necessary to unpin the breather. However, in contrast 
to the results for Hamiltonian systems, once the breather 
starts to move the velocity is unique (independent of A), 
i.e. the MB is a robust attractor of the dynamics. Once 
a induced MB is generated at a value of C, this solution 
can be continued by varying C slowly, with almost no 
variation in velocity. The analysis of Poincare sections of 
this MB shows clearly a quasiperiodic behavior. There- 
fore there is not commensurability between the internal 
frequency ujb and the new frequency associated to the 
velocity w,„b = ^irv^h- This behavior precludes the use 
of fixed point methods (like Newton method) based in 
the periodicity of solutions, to find (numerically) exact 
MB solutions. Also, this quasiperiodicity prevent us to 
extend the Floquet analysis to MB. 



For intermediate couplings 0.60 ^ C '^S 0.72 (shadow 
region in the phase diagram) , the breather portrays ran- 
dom motion. For some time interval, the MB moves 
regularly in one direction, then suddenly remains inmo- 
bile (but quasiperiodic) for a while, changing then its 
motion to the other direction, and so on. A plausi- 
ble scenario is that of the occurrence of a crisis ]3l| ] 
at C « 0.72 which destabilizes the regularly moving 
breathers (of positive and negative velocity, respectively) , 
giving rise to a chaotic attractor consisting of intervals 
(of random length) of approximately regular motion, fol- 
lowed by changes of direction. 

As we have already mentioned, aX C — 0.88 the sta- 
tionary breather solution (quasiperiodic) disappears, and 
only the MB solutions survive. We will refer to this 
breather as spontaneous slow MB. Its velocity is approx- 
imately half of the fast MB and shows a great variation 
as a function of C. Moreover, there is a narrow win- 
dow around C = 0.89 where both the slow and fast MB 
coexist. 

Regular motion MB are also stable against variation of 
parameters other than C. By increasing Fac, the veloc- 
ity of the breather increases, showing a very asymmetric 
profile, as shown in Fig. |l^. As we explain below, the ori- 
gin of this asymmetric shape can be explained in terms 
of forward and backward emission of phonons from the 
moving breather, which suffer a Doppler effect. Above a 
critical value of Fac, which is dependent on C , a shock 
wave is formed and the regular motion becomes again 
diffusive. 



B. Emission of phonons 

An important feature of both pinned quasiperiodic and 
moving breathers is the emission of low amplitude linear 
waves (phonons). In both cases the breather tails clearly 
show a complex quasiperiodic behavior in which many 
frequencies are involved. For the moving breather, these 
tails are also markedly asymmetric, due to the transla- 
tional motion. Other authors have investigated the be- 
havior of Hamiltonian breathers when subject to phonon 
scattering note that in our case it is the breather 
itself the source of phonons. 

In order to investigate the phonon emission we have 
computed the power spectrum of Ui{t) for sites suffi- 
ciently far away from the breather center, as given by 
the expression 



S{lo) 



uj{t)e"^*dt 



(19) 



In all spectra, we can observe peaks corresponding to the 
driving frequency and its odd harmonics, as expected. 
We also observe a broad band spectrum corresponding 
to frequencies in the phonon band, with several reso- 
nant peaks. In both cases we can explain those frequen- 
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cies satisfactorily in terms of emission of phonons by the 
breather. 



Figure 12 shows a time snapshot and the corresponding 
power spectra for a particle in the tails of a quasiperi- 
odic pinned breather. In this case, the resonant peaks 
simply correspond to frequencies which are linear com- 
binations of ujb and o^ncw (the second basic frequency of 
the quasiperiodic breather): 



(20) 



In the case of moving breathers, the peaks also cor- 
respond to the frequencies given above, but shifted by 
the Doppler effect since they are emitted by a moving 
source. Their calculation requires thus a little more care. 
We recall that the propagation of phonons is given by 
the dispersion relations ( |l3|) . The frequency of the emit- 
ting source, in the reference frame of the source itself, is 
LUsrc = nuji, + 2'KmVmh- This second frequency appears 
because the breather is moving over a periodic potential 
with velocity Umb- However, in the reference frame of 
the medium (the lattice) , this frequency will be modified 
according to: 



± 27rg([Jtail)Wmb, 



(21) 



which is the well-known Doppler effect, only that the 
medium is dispersive [^^ . In particular, note how the 
wavevector q of the propagated phonon depends on lu, as 
given by Eqs. (^. Therefore Eqs. (|l|) and ^) have 
to be solved self-consistently for q and Wtaii, and then 
the different peaks of the power spectrum can be worked 
out. The agreement of this calculation with the observed 
frequencies (see Fig. |l^) is excellent. 

Finally we remark that for small lattice sizes, and when 
using periodic boundary conditions, tails in front of and 
behind the breather can overlap. In such cases, solutions 
are similar to the so-called nanopterons of Hamil- 
tonian systems, in which the MB appears to move in a 
"sea" of phonons. 



C. Breather Collisions and Thermal Effects 

Since we are dealing with a system in which pinned 
(periodic and quasiperiodic) and mobile breathers co- 
exist for the same parameter values, it seems natural 
to study their stability against collisions between them. 
Very different kind of events appear depending on the 
breather velocity and the initial conditions (initial dis- 
tance between breathers). The faster the breathers are, 
the more likely they are to destroy each other. The main 
result is that, when the breathers survive to the colli- 
sion, the interaction is mediated by the phonons which 
dress the breather. For large velocity, the tail in front 
of the breather is short (the MB is very asymmetric), so 
the breather cores can overlap and we observe they get 
destroyed. For moderate velocities, the opposite occurs: 



the tails are large, and it seems as if they mediated the 
collision, slowing down the breathers and preventing the 
cores to touch. A clear example of this latter case can 
be seen in the figure |lj, where we observe an "elastic" 
collision in which the breathers approach each other up 
to a distance which is the range of the phonon tail. Note 
that this distance (« 150) is much larger than the core 
breather size (« 10). 

In figure |l5| we show some of the multiple behaviors we 
have observed in simulations: (i) Collision between two 
MB at low Fac and thus low velocity, at C = 0.75. Their 
velocity is low enough to not destroy themselves, but to 
create a new multibreather state; (ii) elastic collision of 
two "slow" breathers at C = 0.89; (Hi) collision which 
gives a mobile two-breather state. The last two collisions 
only differ in the initial conditions. Finally, (iv) breather 
annihilation between a "slow" and "fast" breather at C = 
0.89 where they coexist. 

An interesting phenomenon arising from the simula- 
tions above is the formation of multibreather solution 
both pinned and mobile. These breather "molecules" are 
appear to be linked by "phonon bonds" . Surprisingly, 
these multibreather solutions are more robust against 
changes of parameters. For instance, they can be sub- 
ject to larger Fac without losing regular motion and then 
reaching larger velocities. However, a systematic study 
and rigorous characterization of these configurations will 
be left for further publication. 

Finally, we have incorporated in the equations of mo- 
tion (|l|) a random force £,i{t) with < ^i{t) >= and 
< ^i{t)^jit') >= 2TSi^j6{t - t') in order to simulate the 
Langevin dynamics of our system. For low T {T < 10""*) 
breathers solutions (pinned and mobile) are stable in the 
whole range of coupling C, in the sense that localization 
persists. To be precise, we observe that, if we are in a re- 
gion of parameters where the MB exists, the noise always 
induces motion, which is of course stochastic itself. This 
can be understood again in a scenario of crisis induced 
by the thermal noise. For higher T the thermal excita- 
tion of kink-antikink pairs and other breathers masks the 
original breather. 



VI. DISCUSSION AND CONCLUSIONS 

Dissipative discrete breathers (DB's, for short) are 
generic solutions of forced-damped lattices of non lin- 
ear oscillators. Contrary to their Hamiltonian counter- 
parts, which are severely affected by harmonic resonances 
with the phonon band, the intrinsic localization of energy 
in the dissipative case is not easily destroyed by reso- 
nances due to the efficient damping of the radiation away 
from the localization site. The character of attractor of 
the dissipative DB's allows their numerical continuation 
in parameter space with simpler procedures than those 
needed for the continuation of Hamiltonian DB's, and 
their robustness against all kind of small perturbations 
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(including stochastic ones) ensures their observabihty in 
experimental situations. 

Pinned dissipative DB's which are continued from the 
uncoupled limit experience generically different kinds of 
instabilities by localized modes, namely pitchfork (for- 
ward and backward) and Hopf (supercritical and subcrit- 
ical) bifurcations. Pitchfork bifurcations produce DB's 
with broken mirror symmetry, while Hopf bifurcations 
lead to quasiperiodic DB's. In any case it is remark- 
able that the two basic solutions (one-site and two-site 
periodic breather) are found to be continuable (as sta- 
ble or unstable solutions) up to the continuum limit, 
where continuous translational invariance is restored in 
the model and then a Goldstone mode appears in their 
Floquet spectrum. For no t-too- large forcing, tail analy- 
sis as explained in section [II C explains successfully the 



numerical power spectra in sites away from the breather 
center, as well as DB profiles. 

In certain regions of the parameter space, mobile DB's 
occur as attractors for an open set of initial conditions 
(basin of attraction) in phase space. Their velocity is 
determined by the model parameters, and it is slow com- 
pared with the time scale set by the forcing frequency. 
One class of mobile solutions are connected to the exis- 
tence of depinning modes in the Floquet spectrum of pe- 
riodic DB's, when they are in the vicinity of (symmetry- 
breaking) pitchfork bifurcations. This mechanism allows 
the generation of mobile DB's for values of the coupling 
which are surprisingly low. Another class is related to 
the De-stabilization of quasiperiodic pinned breathers, 
and their velocity is even slower than that of the previ- 
ous class. 

Finally, we tested the robustness of breathers by means 
of collisions and the application of stochastic perturba- 
tions (Langevin noise). It is concluded that breathers 
are quite robust in the sense that the localization per- 
sists; however, in regions of parameters where mobile 
DB's exist, the thermal noise induces random motion in 
the breather. 
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FIG. 2. Inverse of the decay length (in logarithmic scale) 
as a function of the couphng parameter C, for a fixed value 
LJ = 3uJb. We mark two values of C, 0.564 and 0.642, which 
correspond to two rather different values of ^, also used in 
the next figure. The rest of the parameters are cub = 0.27r, 
a = 0.02 and Fac = 0.02. 
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FIG. 1. Wave vector q and inverse of the decay length ^ 
as functions of ui for two different values of the damping, 
a = 0.02 (open circles) and the Hamiltonian case a — 0.0 
(filled ones). The coupling parameter C is in both cases equal 
to 0.75. 
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FIG. 3. Breather profiles (two-site breathers) at a given 
time for two different values of C (0.564 filled circles and 0.642 
open ones, see also the previous figure). The inset shows the 
right hand side tails, where we can observe the existence of a 
phonon in the second case, corresponding to the entrance of 
the third harmonic of wj, in the phonon band. The points are 
connected as a guide to the eye. 
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FIG. 4. Floquet spectrum for the one-site breather near 
C = 0.5296 and closely after the first pitchfork bifurcation. 

All eigenvalues are in a circle of radius cxp{—atb/2) except 
two, one of which crosses the unit circle by +1. The lower fig- 
ure shows the profile of the eigenvector corresponding to this 
unstable eigenvalue (both velocity and position components, 
ej, e,). Note that it is antisymmetric and strongly localized. 
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FIG. 5. Profiles of the breather around the area of the 
pitchfork bifurcation. 




-0.5 0.5 

Re(^l) 



0.2 
0.1 

.cT 

-0.1 
-0.2 [ 



0.2 
0.1 



-0.1 
-0.2 



• 

• 


• 




^ ^ 


■ ^ 





100 200 300 400 500 600 700 800 900 
n 



FIG. 6. Floquet spectrum for the Hopf bifurcation at 
C = 0.871. An eigenvalue and its complex conjugate cross 
the unit circle at a non-zero angle in the complex plane. The 
figure at the bottom shows the corresponding eigenvector. 
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FIG. 7. Two-site quasiperiodic breather. Note that parti- 
cles on both sides of the breather are out of phase. The figure 
below shows the power spectrum of one of the central parti- 
cles. The peaks are linear combinations of the two relevant 
frequencies, ujb and Note that only odd harmonics of 

uJb appear. 
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FIG. 8. Scheme showing all bifurcations found in our model 
up to values of C w 1. The most relevant ones axe probably 
the pitchforks which exchange stability between the one-site 
and the two-site breather (Cpi, Cp2). At C = 0.817 the 

two-site breather has a subcritical Hopf bifurcation, con- 
necting it to a quEisiperiodic breather. At C = 0.88 this 
quasiperiodic solution disappears as it turns into a slow mov- 
ing breather. The grey-shaded areas between C = 0.51 and 
C = 0.96 are those where moving breathers, either fast or 
slow, can be found. 
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FIG. 9. Velocity of MB vs. coupling parameter C. Dashed 
area shows the region in which MB have a diffusive motion. 
It is also showed different regions of "slow" and "fast" MB. 
See the text for details. 
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FIG. 10. The upper figure shows a MB with regular motion 
at C = 0.75. Below, a MB with diffusive motion at C = 0.65. 
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FIG. 12. Time evolution for the tail of quasiperiodic 
breather at C = 0.817 (upper figure) and the correspond- 
ing power spectrum (down figure) whose peaks are given by 
linear combinations of two frequencies uJb and uJnew ■ 
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FIG. 11. Instantaneous profile of a MB at C = 0.75 and 
Fac = 0.045. Note the clear asymmetry of the phonon tails 
in front of and behind the breather. 
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FIG. 13. Time evolution for a particle which is passed by a 
MB (upper). Below, power spectra corresponding to the tails 
before and after the passage of the breather. Note how the 
peaks are shifted by Doppler effect. 




FIG. 14. An elastic collision between two slowly moving 
breathers. It is possible to appreciate how the breathers "see" 
each other through their phonon tails. This collision was ob- 
tained at C = 0.89. 
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FIG. 15. Four scenarios of collisions between breathers. We 
represent the traces by plotting the position of the energy 
maximum, (i) Formation of a pinned "molecule" after a col- 
lision of two MB at C = 0.75 and Fac = 0.011. (ii) "elastic" 
collision at C = 0.89 and Fac = 0.02. (iii) Formation of mo- 
bile "molecule" C = 0.89 and Fac = 0.02. (iv) Annihilation 
of a "slow" and a "fast" breather at C = 0.89 and Fac = 0.02. 
In this last case, the traces after the collision correspond to 
linear radiation (phonons). 
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